CONTINUUM PERCOLATION IN HIGH DIMENSIONS 
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Abstract. Consider a Boolean model S in M. d , where the centers are given by a homo- 
geneous Poisson point process with intensity A and the radii of distinct balls are i.i.d. 
with common distribution v. Some numerical simulations and some heuristic arguments 
suggest that the critical covered volume c c d {v), which is the proportion of space covered 
by £ at critical intensity, may be minimal when v is a Dirac measure. 



1. Introduction and statement of the main results 

Introduction. Consider a homogeneous Poisson point process on M. d . At each point of 
this process, we center a ball with random radius, the radii of distinct balls being i.i.d. 
and independent of the point process. The union £ of these random balls is called a 
Boolean model. This Boolean model only depends on two parameters : the intensity A 
of the point process of centers, and the common distribution v of the radii of the balls. 

Denote by A^(V) the critical intensity for percolation in S. We then consider c c d {y\ the 
proportion of space which is covered by £ when A = \ c d [y)- This quantity is called the 
critical covered volume. Some numerical simulations in low dimension and some heuristic 
arguments in any dimension suggest that the critical covering fraction may be minimal 
when v is a Dirac measure, that is when all the balls have the same radius. 

We show that this is not true in high dimensions. We prove this by studying the critical 
covered volume c d (r~ d fi(dr)) of a given measure /x as d tends to infinity. See the definition 
OH]) below for some remarks about this normalization. When /i is a Dirac measure, this has 
been studied by Penrose in jT3]. Our results yield that whenever /i is not a Dirac measure, 
for all d large enough, c d (r~ d fi(dr)) is strictly smaller than c d (r~ d 5i(dr)) = c d (5i). 

When /i is a Dirac measure, the asymptotic behavior of c d (r~ d 5i(dr)) = c d (5\) is 
given by the associated Galton- Watson process. This is due to the fact that geometrical 
dependencies vanish in high dimensions. This is not true anymore in the general case and 
geometry still plays a role. 

The Boolean model. Let us give here a different - but equivalent - construction of the 
Boolean model. Let v be a finiteS measure on (0, +00). We assume that the mass of v 
is positive. Let d > 2 be an integer, A > be a real number and £ be a Poisson point 
process on M. d x (0, +00) whose intensity measure is the Lebesgue measure on M. d times 
\v. We define a random subset of M. d as follows: 

£(Az/)= |J B(c,r), 

(c,r)6£ 

where B(c, r) is the open Euclidean ball centered at c 6 M d and with radius r e (0, +00). 
The random subset £(Ai/) is a Boolean model driven by \v. 



1. There is no greater generality in considering finite measures instead of probability measures ; this 
is simply more convenient. 
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Critical intensity. We say that E(Az^) percolates if the connected component of E(A^) 
that contains the origin is unbounded with positive probability. This is equivalent de the 
existence of an unbounded connected component of E(Ai/) with probability one. We refer 
to the book by Meester and Roy [TT] for background on continuum percolation. Set: 

X c d {v) = inf{A > : £(Ai>) percolates}. 

One easily checks that A^(z^) is finite. In [6] it is proven that \ c d {v) is positive if and only if 



(1) 



J r d v[dr) < +oo. 



We assume that this assumption is fulfilled. 

Critical intensity as a function of u, critical covered volume. By ergodicity, the 
Boolean model E(Af) has a deterministic natural density. This is also the probability 
that a given point belongs to the Boolean model and this is given by : 

P(0 e S(Ai/)) = 1 - exp ^-A J v d r d v{dr) 

where v d denotes the volume of the unit ball in IR d . The critical covered volume c d [y) is 
the density of the Boolean model when the intensity is critical : 



cS(f ) = 1 - exp \-\ c d {v) J v d r d v(dr 

It is more convenient to study the critical covered volume through the normalised thresh- 
old defined by : 



K(y) = K(y) J v d (2r)*v(dr). 

Note : 

«-) = i-«p(-^ 

The factor 2 d may seem arbitrary here. Its interest will appear in the statement of the 
next theorems. 

We will now give two scaling relations. For all a > 0, define H a v as the image of v 
under the map defined by x h- > ax. By scaling, we get: 

(2) X d (H a u) = X d (u). 

This is a consequence of Proposition 2.11 in [TT], and it may become more obvious when 
considering the two following facts : a critical Boolean model remains critical when 
rescaling and the density is invariant by rescaling ; therefore the critical covered volume 
and then the normalised threshold are invariant. One also easily checks the following 
invariance: 

(3) X d (au) = X d (u). 

That is because of these two invariance relations that we prefer to study the normalized 
critical intensity X d (v) instead of the critical intensity X c d (v) 

It has been conjectured by Kertesz and Vicsek [9] that the critical covered volume - or 
equivalently the normalized critical intensity - should be independent of u, as soon as the 
support of v is bounded. Phani and Dhar [T] gave a heuristic argument suggesting that 
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the conjecture were false. A rigorous proof was then given by Meester, Roy and Sarkar 
in [12]. More precisely, they gave examples of measures v with two atoms such that: 



(4) A5(i/) > AS(*i). 

As a consequence of Theorem 1.1 in the paper by Menshikov, Popov and Vachkovskaia 
[13], we even get that \%{y) can be arbitrarily large. More precisely, if 

71-1 

(5) v{n,a) = ^a dk 5 a - k . 

k=0 

then 

(6) X d (u(n,a)) — > nX d (6i) as a — > oo. 

Actually the result of [13] is the following much stronger statement: A^(z/(+oo, a)) — > 
X c (Si) when a — > oo. The convergence ([6]) is implicit in the work of Meester, Roy and 
Sarkar in [12], at least when n = 2. There were also heuristics for such a result in [I]. 

By Theorem 2.1 in [6], we get the existence of a positive constant Cd, that depends 
only on the dimension d, such that: 

ASM > c d . 

To sum up, (•) is not bounded from above but is bounded from below by a positive 
constant. In other words, the critical covering volume c c d {-) G (0, 1) can be arbitrarily 
close to 1 but is bounded from below by a positive constant. Our investigations were 
motivated by the following question: 

(7) Is the critical covered volume c c d (v) minimized when v — Si ? 

In the physical literature, it is strongly believed that the answer is positive at least 
when d = 2 and d = 3. This conjecture is supported by numerical evidence (to the best 
of our knowledge, the most accurate estimations are given in a paper by Quintanilla and 
Ziff [15] when d = 2 and in a paper by Consiglio, Baker, Paul and Stanley [2] when 
d — 3) . On Figure [H we plot the critical covered volume in dimension 2 as a function of 
a and for different values of p when v — (1 — a)b\ + ap~ 2 5 p . The data for finite values 
of p come from numerical estimations in [15] , while the data for infinite p come from the 
study on the multi-scale Boolean model. See section 1.4 in [8] for further references. It 
is also supported by some heuristic arguments in any dimension (see for example Dhar 
[3]). See also [1]. In the above cited paper [12], it is noted that the rigorous proof of (jl]) 
suggests that the answer might be positive for any d > 2. 

In this paper, we show that the answer is actually negative for all d large enough. 

Critical intensity in high dimensions : the case of a deterministic radius. 

Assume here that the measure v is a Dirac mass at 1, that is that the radii of the balls 
are all equal to 1. Penrose proved the following result in [T3] : 

Theorem 1.1 (Penrose), lim A^(<5i) = 1. 

d—^oo 

With the invariance properties of A^, this limit can readily be generalized to any con- 
stant radius : for any a, a > 0, 

lim X d (a5 a ) = lim \ c d {5i) = 1. 

d—too d—too 
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Figure 1. Critical covered volume as a function of a for different values 
of p. From bottom to top: p = 2, p = 5, p = 10 and the limit as p — > oo. 



Theorem 11.11 is the continuum analogue of a result of Kesten [10] for Bernoulli bond 
percolation on the nearest-neighbor integer lattice Z d , which says that the critical perco- 
lation parameter is asymptotically equivalent to l/(2d). 

Let us say a word about the ideas of the proof of Theorem 11.11 

The inequality A^(5i) > 1 holds for any d > 2. The proof is simple, and here is the 
idea. We consider the following natural genealogy. The deterministic ball -8(0, 1) is said 
to be the ball of generation 0. The random balls of £(A5i) that touches -8(0, 1) are then 
the balls of generation 1. The random balls that touches one ball of generation 1 without 
being one of them are then the balls of generation 2 and so on. Let us denote by the 
number of all balls that are descendants of -8(0, 1). There is no percolation if and only if 
Nd is almost surely finite. 

Now denote by m the Poisson distribution with mean Xvd2 d '■ this is the law of the 
number of balls of E(A#i) that a touch a given ball of radius 1. Therefore, if there were 
no interference between children of different balls, Nd would be equal to Z, the total 
population in a Galton- Watson process with offspring distribution m. Because of the 
interferences due to the fact that the Boolean model lives in M. d , this is not true : in fact, 
Nd is only stochastically dominated by Z. Therefore, if \vd2 d < 1, then Z is finite almost 
surely, then is finite almost surely and therefore there is no percolation. This implies 

AS(tfi) = v d 2%(5 1 ) > 1. 
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The difficult part of Theorem 11.11 is to prove that if d is large, then the interferences 
are small, then N d is close to Z and therefore there is percolation as soon as X d (5i) = 
v d 2 d X c d (S) > I. 

To sum up, at first order, the asymptotic behavior of the critical intensity of the Boolean 
model with constant radius is given by the threshold of the associated Galton- Watson 
process, as in the case of Bernoulli percolation on Z d : roughly speaking, as the dimension 
increases, the geometrical constraints of the finite dimension space decrease and at the 
limit, we recover the non-geometrical case of the corresponding Galton- Watson process. 

Critical intensity in high dimensions : the case of true random radii. If /x is a 

finite measure on (0, +00) and if d > 2 is an integer, we define a measure [i d on (0, +00) 
by setting : 

(8) fi d (dr) = r- d fi(dr). 

Note that, for any d, the assumption ([1]) is fulfilled by fj, d , and that (Si)d = S\. Note also 
that \id is not necessarily a finite measure. However the definitions made above still make 
sense in this case and we still have X d c {\Xd) £ (0, +00) thanks to Theorem 1.1 in [7]. 

We will study the behavior of \ d {fid) as d tends to infinity. Let us motivate the definition 
of Hd with the following two related properties : 

(1) Consider the Boolean model E(A/x ( j) on M. d driven by Xfi d where A > 0. For any 
< s < t < 00, the number of balls of T,(Xfj, d ) with radius in [s,t] that contains 
a given point is a Poisson random variable with intensity: 

/ v d r d Xjj,d(dr) = v d X/i([s : t\). 

J[s,t] 

Loosely speaking, this means that contrary to what happens in the Boolean model 
driven by A/x, the relative importance of radii of different sizes does not depend 
on the dimension d in the Boolean model driven by A/x^. 

(2) A closely related property is the following one. Consider for example the case 
/x = a5 a + (35},. Then, a way to build E(A/Xd) is to proceed as follows. Consider 
two independent Boolean model: T, A , driven by XaSi, and S B , driven by Xf35i. 
Then set £(Axx d ) = aS A + 6S B . 

We prove the following result : 

Theorem 1.2. Let /x be a finite measure on (0,+oo). We assume that the mass of /x is 
positive and that /x is not concentrated on a singleton. Then : 

limsup - In (x d (fi d )) < 0. 

Note that (Si) d = S±. A straightforward consequence of Theorem 11.21 and Theorem 
11.11 - or, actually, of the much weaker and easier convergence of ln(A^(5i)) to - is the 
following result: 

Corollary 1.3. Let /x be a finite mesure on (0, +00). We assume that the mass of /x is 
positive and that /x is not concentrated on a singleton. Then, for any d large enough, we 
have: 

X d (Hd) < ArfOi), or equivalently c c d (fi d ) < c c d {5i). 

In fact, Theorem 11.21 follows from the particular case of radii taking only two different 
values, for which we have more precise results. 

5 



Critical intensity in high dimensions : the case when the radii can take two 
values. To state the result, we need some further notations. Fix p > 1. Fix k > 1. 

Set r x = r k+l = 1 + p, and for i G {2, ...,&}, r* = 2. For (a^Kfc+i G [0, l) fc , we 
build an increasing sequence of distances (di)i<i<fc+i by setting di = 1 + p and, for every 
ie{2,...,k + l}: 

d\ = d\_ x + 2r i a i d i _ 1 + r t 2 . 

Note that the sequence (dj)i<i<fe+i depends on p, fc, and the a^s. 
We denote by P(a2, . . . , a^+i) = dfc+i- Now set, for every k > 1, 



(9) Kp(^) = inf max 



4p \ 2p 

/ 




Finally, let: 

(10) ^ = inf^). 



We give some intuition on k£ in Section I2~T| make a few further comments and simulations 
for K c p (l), K c p (2) and ^(3) in Section E2J We will see that K c p G (0, 1) for all p > 1. This 

^ gives the asymptotic behavior of \ c d (pd) when p charges two distinct points : 

Theorem 1.4. Let b > a > 0, a > and > 0. 5et p = a5 a + anc? p = b/a > 1. 
Then 

(11) Jjmiln(AS W )=ln(^)<0. 

Remember that in the case of a deterministic radius (assumptions of Theorem II. lft . the 
first order of the asymptotic behavior of the critical intensity in high dimensions is given 
by the threshold of the associated Galton- Watson process. 

In the case of two distinct radii (assumptions of Theorem ll.4p . the associated Galton- 
Watson process is two-type, one for each radius. At first order and on a logarthmic scale, 
the asymptotic behavior of the critical intensity in high dimensions is given by k c , and 
we could expect that, as in the deterministic would be given by the threshold of 

this two-type Galton- Watson process. We will see in Lemma 12.11 that it is actually the 
case, but only when p is not too far from 1. When p is large, this is no longer the case 
and geometry still plays a role. See the discussion in Section 12.21 

The proof of Theorem 11.41 is given in Section El The main ideas of the proofs are give 
in Section I2TT1 Theorem 11.21 is an easy consequence of Theorem II .41 The proof is given in 
Section [3j Corollary 11.31 is a straightforward consequence of Theorems 11.11 and 11.21 We 
give some intuition and results about n c p is Section 12.11 and 12.21 

2. The case when the radii take two values 

2.1. Notations and ideas of the proof. In the whole proof, we fix p > 1 and k > 0. 

Once the dimension d > 1 is given, we consider two independent stationary Poisson 
point processes on M d : Xi an d Xpi with respective intensities 

Ai = 7 and A, 



v d 2 d p v d 2 d p d ' 

With xi an d Xpi we respectively associate the two Boolean models 

E x = (J B(x, 1) and S p = (J B(x, p). 



Note that E p is an independent copy of pEi. Note also that the expected number of balls 
of Ei that touches a given ball of radius 1 is n d . Thus the expected number of balls of 
E p that touches a given ball of radius p is also n d . 

We focus on the percolation properties of the following two-type Boolean model 

E = Ei U E p . 

We begin by studying the existence of infinite /c-alternating paths. For k > 1, an infinite 
^-alternating path is an infinite path made of balls such that the radius of the first ball 
is p, the radius of the next k balls is 1, the radius of the next ball is p and so on. For 
a fixed k > 1, we wonder whether infinite /c-alternating paths exist and seek the critical 
threshold K p {k) for their existence. A natural first step is to study the following quantities 

(12) N = #{xi€ X p: INI < 2p}, 

{%k+\ £ Xp '■ ^( x i)i<i<k G Xi distinct such that 
||xi|| < 1 + p, Vi £ {1, . . . , k — 1} \\xi + i — Xi\\ < 2, 
\\x k+ i - x k \\ < 1 + p 

Fix k > 1. Remember that n p {k) is defined in (Q. 

A lower bound for K p (k). In Subsection 12. 3[ we obtain lower bounds for K c p (k) by looking 
for upper bounds for E(N k ). On one side, a natural genealogy is associated to the 
definition of N k (see also the comments below Theorem 1 1 . 1 1 and in Section [2T2T) . We start 
with an ancestor xq located at the origin. We then seek his children in xi H B(xo, 1 + p) 
: they constitute the first generation. If x\ is one of those children, we then seek the 
children of x\ in xi H B(x\,2) to build the second generation and so on. On the other 
side, the process lives in and the geometry induces dependences: if x\ and x[ are 
two individuals of the first generation, their children are a priori dependent. If we forget 
geometry and only consider genealogy, we get the following upper bound: 

E{N k )<X 1 \B{- i l + p)\ (n Al l 5 (-' 2 )l) X p \B(-,l + p)\. 

But the points of the last generation are in B(0, 2p + 2k). So if we forget genealogy and 
only consider geometry we get the following upper bound: 

E{N k ) < \ p \B(0,2p+2k)\. 

Expliciting the two previous bounds and combining them together, we get: 

E(N k ) < max K - ^, ; 

V 4p p 

In this upper bound, the first argument of the maximum is due to genealogy while the 
second one is due to geometry. To get the geometrical term, we considered the worst case: 
the one in which, at each generation 2, X{ IS cLS far from the origin as possible This gives a 
very poor bound. To get a better bound, we proceed as follows. Fix a 2 , . . . , a k+ \ G [0, 1). 
As before, we set T\ = r k+1 = 1 + p, and for i £ {2, ...,&}, T{ = 2 and we build 
the increasing sequence of distances (rfi)i<i<fc+i by setting di = 1 + p and, for every 
i £ {2,...,fc + 1}: 

d l. = + 2r i a i di- 1 + r\. 
See Figure |2] for a better understanding of these distances d\. 
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Figure 2. Definition of the distances d{. Circles in plain line are of radius 
1 + p and 2. 



Denote by N k (a 2 , ■ ■ ■ , Ofc+i) the number of points x k+ i E Xp f° r which there exists a 
path xi, . . . ,x k fulfilling the same requirement as for N k and such that ~ di for all 
i. Proceeding as before, we obtain the following upper bound: 



d 

fc+ i (1 + pf / rr M ^ KVfa, • • • , afc+i) > 
E(N k (a 2 , • • • , a k +i)) < mm « + 




Here again, the first argument of the maximum is due to genealogy while the second one 
is due to geometry Optimizing then on the a^'s, we get: 

t?(at \<t ■ I fc+i( 1+ P) 2 / TT ft ^ K%,...,a w 
E(N k )< sup mm I k — / 1 - a, 2 , 



a 2 ,...,a k+1 



A precise statement is given in Lemma 12751 The precise value of the threshold Kp(k) given 
in (Q is then the value such that the above upper bound converges to when k < n c p {k). 
This heuristic will be precised in Subsection 12.31 we will prove there that when k < n c p {k), 
E(N k ) converges to as d tends to infinity, and this will imply that there exists no infinite 
/c-alternating path. 
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An upper bound for K p (k). If, on the contrary, k > n p {k), then we will prove that E(N^) 
does not converge to 0. Actually, to prove that when k > K p (k), there exist infinite k- 
alternating path, we will show, in Subsection \2A\ the following stronger property : With 
a probability that converges to 1 as d tends to infinity, we can find a path which fulfills 
the requirements of the definition of N/. - or more precisely of Nk(a,2, ■ ■ ■ , cik+i) for some 
02, ... , czfc+i nearly optimal - and which fulfills some extra conditions on the positions of 
the balls. This is Proposition 12.91 and this is the main technical part of this paper. Those 
extra conditions provide independence properties between the existence of different paths 
of the same kind. We can then show the existence of many such paths and concatenate 
some of them to build an infinite fc-alternating path. Technically, the last step is achieved 
by comparing our model with a supercritical oriented percolation process on Z 2 . In this 
comparison, an open bond in the oriented percolation process corresponds to one of the 
above paths in our model. This comparison with oriented percolation was already the 
last step in the paper of Penrose [H] • 



From infinite k-alternating paths to infinite paths. Recall k c = inffc>i K c {k). With the 
previous results, it rather easy to show that there is no percolation for d large enough as 
soon as k < k c . When k > n c p then k > K c {k) for a k > 1. Therefore there is /c-alternating 
percolation and therefore there is percolation. 

2.2. The threshold K p and the associated Galton- Watson process. We begin with 
a couple of results about K p (l) and K p . 

Lemma 2.1. Let p > 1. Then 
• < K p (l) < 1. More precisely : 

Ifl<p<2 then K °Jl) = while if p > 2 then K c Jl) ~ ^ 



1 + p 7 J r ~ pw 1 + p 



0<K c p <l. 



There exists po > 2 such that if p < po, then k c = k c (1). This implies 



If 1 < p < 2 then K p = - v ^ , while if 2 < p < p then K c p = — — 



1 + J p 1 + p 

On Figure[3l we plot n p (i), for i G {1,2, 3}. The data come from numerical estimations. 



Proof. • By definition, n c (l) = inf max(0i(a), 02(o)) 5 where 0i,02 : [0,1) — > K are 

0<a<l 

defined by: 

, , , 2^p P V2 
M a ) = 7^, — rr, ^ttt and Ma) ~ 



(l + p)(l-a2)V4 ™ , (1 + p)v ^ T7i - 

If p < 2 then 0i(O) > 02(0). As 0i is increasing and 02 is decreasing, we get: 

K p (l) = inf Ma) = 0i (0) = 

^ 0<a<l 1 + p 

Assume, on the contrary, p>2. Set 

p 2 -4 



p 2 + 4 
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€ [0,1). 




Figure 3. Simulations for K p (i), iG {1,2,3}, from left to right. 
Then 0i(a) = 02(a)- As <pi is increasing and cj>2 is decreasing, we get: 



J A + d 1 

K c p (l) = Ma) = Ma) ~ 



l + P 



Clearly we have, for every k > 1 : 
Kp(k) > _ inf 



i/(fc+i) 

4p \ / 4p 



°- a2 -" afc+i<1 V(i+p)Vn2<i< fc+ i(i -«?))/ V(1+p) ' 

Therefore, as = inf fc >x n c p (k), 



i/(fc+i) 



0< ((TT^) I/2 ^^^ (1)<h 

The last inequalities imply that = k c (1) if 1 < p < 2. In fact, as soon as 

4p ^ 1/3 



(13) 

it is true that k c (1) < n c p (k) for all k > 2 and therefore that k° — k c (1). As the inequality 
in f fT3|) is strict for p = 2, we obtain by continuity the existence of p > 2 such that for 
every p G (1, p ),n c p = n c p (l). □ 

The associated Galton- Watson process. The associated Galton- Watson process is a two 
type Galton- Watson process, one type for each radius. Consider for example the offspring 
distribution of type p of an individual of type 1. We define it to be the number of balls of 
a Boolean model directed by X P S P that intersects a given ball of radius 1. Therefore, this 
is a Poisson random variable with mean X p Vd(l + p) d . The other offspring distribution 
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are defined similarly. We moreover assume the the offspring of type 1 and p of a given 
individual are independent. 

Consider the matrix of means of offspring distributions: 

_ (\ lVd {l + l) d \v d (l + P) d 
d ~ \\xv d (l + p) d \Mp + pY 

For example, the top-right term is the mean offspring of type p of an individual of type 
1. We have: 

i (V* 

,<2 



M d = n d I 



Let r d denotes the largest eigenvalue of M d . The extinction probability of the Galton- 
Watson process is 1 if and only if r d < 1. We have: 

<i + pY c 



If we fix k and let d tends to infinity, the critical k for the Galton- Watson process is thus 

K = . 

1 + p 

With Theorem 11.41 and Lemma 12.11 we thus see that the comparison with the Galton- 
Watson is asymptotically sharp on a logarithmic scale only when 1 < p < 2. 

2.3. Subcritical phase. Let p > 1 be fixed. We consider, in M. d , the two-type Boolean 
model E introduced in Subsection I2.1[ with radii 1 and p and respective intensities 

~d ,-d 
Vd2 d P = Vd2 d p d 



Ai = — —7 and A, 



depending on some k G (0,1). The aim of this subsection is to prove the following 
proposition: 

Proposition 2.2. Let p > 1 be fixed. If k < K c p , then, as soon as the dimension d is large 
enough, percolation does not occur in the two-type Boolean model E. 

In the following of this subsection, we fix p > 1 and < k < K c p . 

We start with an elementary upper bound, in which we do not take into account the 
geometrical constraints. We recall that the N' k s have been introduced in (I12|) . 



Lemma 2.3. £(iVo) = n d and, for k > 1, E(N k ) < 



K fc+1 (l + p) 2N 
4p 



Proof. The result for iVo follows directly from the equality E(Nq) = X p \B(0,2p)\. 
Take now k > 1. We have : 

(14) E{N k ) < \ X \B(; 1 + p)\ \ X \B{-, 2)| j \ p \B{-, 1 + p)\ 
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where B(-,r) stands for a ball with radius r and center unspecified. This can for instance 
be seen as follows : 

E(N k ) < E I ^2 l Xl eB(o,i+p) ■ ■ ■ lx k+1 eB(x h ,i+ P ) 

\xi,...,x k £xi distinct, x k+1 £Xp 
= ^l^p / dX\ . . . dXk+ll Xl £B{0,l+p) ■ ■ ■ ^-x k+1 eB(x k ,l+p), 

jRd(k+l) 

which gives (fl4|) . The lemma follows. □ 

To give a more accurate upper bound for the AVs, we are going to cut the balls 
into slices and to estimate which slices give the main contribution. For x G M. d \ {0}, 
< a < b < 1 and r > 0, we now define : 

If a > : B(x, r, a, b) = |y G M. d : \\y — x\\ < r and ar < — x, | — < br j- , 

If a = : B(x, r, a, b) = |y G M d : \\y — x\\ < r and (^y — x, — br 

The next lemma gives asymptotics for the volume of these sets: 
Lemma 2.4. For x G M d \ {0}, < a < b < 1 and r > 0, 



1 n „ ( \B(x, r, a, b)\ 

Actually we will only use 



lim — In ( - — v ' ' — ) = In ( ryl — a? 



lim sup — In ( - — — ' Q ' ^ | < l n fry/T 
d^+oo d \ v d ) V 



a 2 



Proof of Lemma 12.41 Note that it is sufficient to prove the lemma for x = ei, first 

vector of the canonical basis, and r = 1. 

First, if a = 0, the result follows directly from the inequality v d /2 < \B(ei, 1, 0, b)\ < v d . 
Assume next that a > 0. On the one hand, B(ei, 1, a, b) is included in the cylinder 



{x = Oi, . . . , x d ) : Xi G [a, 1] and || (0, x 2 , ■ ■ ■ , x d ) \\ < Vl - a 2 }, 
which implies 



(15) ^(d, 1, a, 6)| <v d -Wl-a 2 (1-a). 

On the other end, by convexity, B(ei,l,a,b) contains the following difference between 
two homothetical cones: 

\ — 

{x = Oi, • • • , Xd) : X\ G [a, b] and ||(0,x 2 , ...,x d )\\ < Vl 7 ^— — }, 
which implies 

(16) Vd ~ lVT d^ K 1 -«)-(!- b )) < l 5 K 1. 6)1- 

The lemma follows from flT5l) and (1161). □ 



We can now improve the control given in Lemma 12.31 
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Lemma 2.5. For every k > 1, 

lim sup -hx(E(N k )) 

< In sup nun I k — / 1 - aj), 

\o<a 2 ,..., ak+1 <i y 4 P y 2 <y + i 2 p 

Proof. • Fix N > 1. Note that the ball B(x,r) is the disjoint union of the slices 
B(x,r,n/N, (n + 1)/N) for n 6 {0, . . . , N - 1}. For any n 2 , . . . ,n k+1 G {0, . . . , N — 1}, 
we set 

nj , , + 1 
at = — and aT = — ^— • 

We focus on the contribution from a specific product of slices : 

{x k +i G x P '■ 3(xj)i<i< fc G xi distinct with 
||xi|| <l + p, Vz G {1, . . . , k - l}x i+ i G B(xi,2,a i+1 ,af +1 ) 
x k+1 G B(x fe , 1 + p, ajfc+i, a£ +1 ) 

Then we have : 

(17) N k < ^A fc (n 2 ,...,n fc+1 ), 

where the sum is over (n 2 , . . . , n k+1 ) G {0, . . . , N — l} k . 

• As we can check that the points contributing to N k (n 2 , ■ ■ ■ , n k+ \) are in B(0, 1>(a\ , . . . , a k+1 )), 
we get : 

E(N k (n 2 , . . . , n k+1 )) < \ p v d V(a%, a^ +1 ) d , 

this leads to : 

(18) limsup \ In {E{N k {n 2 , n k+1 ))) < In f KV ( a t,---> a k+i, 



• Besides, proceeding as in the proof of Lemma 12. 3[ we obtain : 
E(N k (n 2 ,...,n k+1 )) < X 1 \B(0,l + p)\ ( JJ Ai|5(-, 2, a h a+)\ J A p |£(-, 1 + p, a fe+1 



,i=2 



With Lemma T2 .41 we deduce : 



(19) limsup i In E(iV fc (n 2 ,...,n fc+ i)) < In l^^-tftj J] (1 



a!— >oo 



2<j<fc+l 



• From (ED, (J32D and (p2D we finally get : 
lim sup 

d— >+oo u> 

< ln( max min ( ^S±fl I TT (1 _ 2) ^(4, ■ ■ ■ , 

- \« -..in v » ^ 4p y 2 ^| +1 v 2p 

As T> is uniformly continuous on [0, l] fe , we end the proof by taking the limit when N 
goes to +oo. □ 

The next step consists in taking into account all k > simultaneously; we thus intro- 
duce 
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(20) 

(( y G x P '■ 3k > 1, 3(x i )i< i < fc G xi distinct with 

^ \\xi\\<l + p, Vie{l,...,A;-l}||x i+ i-x < || <2, )> U G x P : IMI < 2p} 

[ jjy - x k \\ < 1 + p 

Lemma 2.6. 7/k < K c p , then limsup - \n(E(N)) < 0. 
Proof. We have : 

(21) E(N)<Y,E(Nk)- 

k>0 

As k < Kp, Lemma [2.51 ensures that for every k > 1 : 

(22) limsup - ln(8(iV fc )) < 0. 

d—>+oo Q> 



Moreover, the assumption k < n c p also implies, thanks to Lemma I2TT] that k < 1. We can 
then choose /co large enough to have : 



Ap 



< exp(-l). 



With Lemma |2.3[ we thus get : 

E(N ) + E ( N k) < K<i + exp(-d) kM = Kd + exp(-d) ^3^- 

fc>fc fe>0 

With (J2TJ) and fl22]), this ends the proof. □ 

The next lemma is elementary 
Lemma 2.7. Assume k < 1. 

T/ien i/ie connected components of [^J -B(x, 1) are bounded with probability 1. 

Proof. For any integer A; > 0, denote by M(k) the number of balls with radius 1 linked 
to -8(0, 1) by a chain of k distinct balls with radius 1. Proceeding as in the proof of 
Lemma I2.3[ we get : 

E(M k ) < K di - k+l \ 

Now denote by M the number of balls with radius 1 linked to -8(0, 1) by a chain of 
(perhaps no) balls with radius 1. Then : 

E(M)<E[J2M h )^ I ^- d <+oc. 

\k>0 / 

Therefore, M is finite with probability 1. So the connected components that touch -8(0, 1) 
are bounded with probability 1. So with probability 1, every connected component is 
bounded. □ 

Proof of Proposition I2T2T . Remember that we proved in Lemma 12.11 that k c < 1. 
Take k such that < k < k c < 1. 

Let £i be the set of random balls with radius p that can be connected to -8(0, p) through 
a chain of random balls with radius 1 (we consider the condition as fulfilled if the ball 
touches -8(0, p) directly). Let £2 be the set of random balls with radius p that are not in 
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£1, but that can be connected to 5(0, p) through a path of random balls in which there 
is only one ball with radius p. We define similarly £3,^4 and so on and denote by £ the 
disjoint union of all these sets. 

We have #£1 = N. (Remember that N has been defined in (120]) .) By Lemma [2.6[ we 
have : 

limsupiln(^i)) < 0. 

d— !-+oo Oi 

Take some p > and assume from now on that d is large enough to have 

i]n(tf(#fc))<-/*. 

Pour tout k > 1, on a alors : E(#£ k ) < (E^^)) 1 * < exp(-dkp). 

As £ = Ufc>i£fc, we deduce from the previous inequalities that £ is finite with probability 
1. So any unbounded connected component of Si U E p that touches -8(0, p) is included in 
Si. As k < 1, Lemma [2~71 rules out the possibility of an unbounded connected component 
in Si. So with probability 1, the connected components of Si U S p that touch B(0,p) 
are bounded, which ends the proof. □ 

2.4. Supercritical phase. We fix here p > 1. We consider once again the two-type 
Boolean model S introduced in Subsection 12. II and we fix an integer k > 1. 

For every n > 0, we set R n = p if k + 1 divides n and R n = 1 otherwise. We say that 
percolation by fc-alternation occurs if there exists an infinite sequence of distinct points 
(a: n ) ne N in M. d such that, for every n > 0: 

- x n E XR n - 

- B(x n , R n ) n B(x n+1 , R n+1 ) ^ 0. 

In other words, percolation by /c-alternation occurs if there exists an infinite path along 
which k balls of radius 1 alternate with one ball of radius p, ie if there exists an infinite 
/c-alternating path. The aim of this subsection is to prove the following proposition: 

Proposition 2.8. Let p > 1 and k > 1 be fixed. Assume that n G (K c p (k),l). If the 
dimension d is large enough, then percolation by k-alternation occurs with probability 
one. 

As announced in subsection 12.11 percolation by /c-alternation of the two-type Boolean 
model in the supercritical case will be proved by embedding in the model a supercritical 
2-dimensional oriented percolation process. 

We thus specify the two first coordinates, and introduce the following notations. When 
d > 3, for any x G lR d , we write 

x = (x',x") G M 2 x R d ~ 2 . 

We write B'(c, r) for the open Euclidean balls of M 2 with center c G M 2 and radius r > 0. 
In the same way we denote by B"(c,r) the open Euclidean balls of IR d-2 with center 
c G IR d ~ 2 and radius r > 0. 

2.4.1. One step in the 2-dimensional oriented percolation model. The point here is to 
define the event that will govern the opening of the edges in the 2-dimensional oriented 
percolation process : it is naturally linked to the existence of a finite path composed of k 
balls of radius 1 and a ball of radius p, whose positions of centers are specified. 
We define, for a given dimension d, the two following subsets of M. d : 

W = cT 1/2 ((-1,1) x (-1,0) x R d - 2 ) , 

W + = d~ 1/2 ((0,0) x (0,1) x R d ~ 2 ) . 
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For x E W we set : 

(23) Q + {x) — i There ex ist distinct £i, . . . , Xk E Xi H ^ + an d x k+i E Xp H 1 
\ such that Xo, Xi, . . . , x^+i is a path J ' 

Our goal here is to prove that the probability of occurrence of this event is asymptotically 
large : 

Proposition 2.9. Let p > 1 and k > 1 be a fixed. Assume that k E (i% c p (k), 1) and choose 
p E (0, 1). // the dimension d is large enough, then for every x E W , 

P(Q + (x)) > P . 

Note already that by translation invariance, P(Q + (x)) does not depend on x", so we 
can assume without loss of generality that x" = 0. In the sequel of this subsection, p > 1 
and k > 1 are fixed. 

We first recall the definitions of the (di)i<i<fc+i and of K c p (k) we give in the introduction. 
We set ri = Tfc + i = 1 + p and for i E {2, . . . , k}, r, = 2. Then, for a given sequence 
(otj)2<i<fc+i 6 [0, l) fc , we build an increasing sequence (di)i<i<jt+i by setting di = 1 + p 
and, for every i E {2, . . . , k + 1} : 

Finally, we note £>(a 2 , . . . , Qfe+i) = and we set 

4p \ 2p 



= m ^ max 

0<a2,...,a fe+1 <l 



^ V 1 + -°) 2 ^ruw 1 - a ?) / ' PK • • • ' ak+l> ) 



The first step consists in choosing a nearly optimal sequence (aj)2<j<fc+i E [0, l) k satisfying 
some extra inequalities : 



fc+1 



Lemma 2.10. We can choose (aj)2<i<fc+i E [0, 1) such that 

y 2<j<fc+i r 

Proof. As Kp(k) < k, we can choose (a°)2<i<fc+i G (0, l) k such that the two following 
conditions 

i 

fc+i 

4p 



2p 



(25) k > 

(26) k > 

dk+i 

are fullfilled for (aj)j = (a°)j. We fix (a 3 , • • • , afc+i) = (a°, • • • , a° +1 ). Note that 

jf ". 0/2 ' — ^ k fc+1 is continuous and increasing, 
2p 



g : a 2 ^ n k+1 - — "|~ / J^J (1 — a?) is continuous and decreasing, 

P Y 2<j<fc+l 

and that lim a2 _>.i #(02) = 0. Moreover, Conditions f[26]) and (1251) ensure that /(a°) > 1 
and ^(a°) > 1. Thus we only require that /(a^) > fi'( a 2) to prove the lemma. 
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If /(a^) — g{ a V)i we can take a 2 > a 2 such that 1 < g(a 2 ) < /(a^) : then 7(02) — 
f( a< 2) > 9{ a 2) > 1 and the lemma is proved. □ 

Note that((2lD implies (125|) and (|26|). 

As explained in subsection 12. 3[ the main contribution to the number Nf. of centers Xk+i 
of balls of radius p that are linked to a ball of radius p centered at the origin by a chain 
(xi)i<i<k of k balls of radius 1 - see the precise definition (fl2|) - is obtained for \\xi\\ ~ di, 
where the d[s are build from a (nearly) optimal sequence {ai)2<i<k+i £ [0, So we fix a 
nearly optimal family (ai)2<i<k+i £ (0, l) fc satisfying f[24]) . we build the associated family 
of distances (d;)i<j<fc+i G (0, l) fc and we are going to look for a good sequence of centers 
(^)i<i<fc+i e (0, l) fc with \\xi\\ ~ di. 

We thus introduce the following subsets of M. 2 : 

D = (-cT 1 / 2 ,^ 1 / 2 ) x (-d-V^o), 
Vi<={l,...,A; + l} D; = (O,^ 1 / 2 ) 2 , 
and the followining sets in M d ~ 2 

Co = {0}, 

Vie{l,...,fc+1} G'l = B"(0,d i -2d~ 1 )\B"(0,d i -3d~ 1 ). 

Finally, for % G {0, . . . , k + 1}, we set Ci = D\ x C". Note that for d large enough, these 
sets are disjoint. The next lemma controls the asymptotics in the dimension d of the 
volume of these sets 

Lemma 2.11. For every i G {1, . . . , k + 1} : 

lim — In - — — = lim — In - — — = In di . 

d->+oo d Vd-2 d^+oo d Vd 

Proof. This can be proven by elementary computations. □ 

Each Xi will be taken in Ci, but we also have to ensure that the (x»)i<i<jb+i form a 
path. Note that for i G {2, . . . ,k + 1}, we have < + a^i < di, which legitimates 
the following definition. See also Figure EJ For % G {2, . . . , k + 1} and d large enough, we 
denote by di the unique real number in (0, 7r/2) such that 

di-i + (XiTi 7—1/2 
COS^j = : h a ' . 

di 

We introduce next, for y G Cj_i, the following subset of M. d ~ 2 : 

D'!{y") = {z" G C>! : (z",y") > \\y"\\.\W% cos ft} 

We also set _D ' = C ' and D"(y") = G'{ for every y £ C . Finally, we define for every 
i G {1, . . . , k + 1} and ?/ G Cj_i : 

A(y) = Di x L>f(n") c a, 

and D = D' x 

Lemma 2.12. • // the dimension d is large enough, for every i G {1, . . . , k + 1} and 

y e Ci_i, 

A(y) c s^^na 

• Let xo G -Do- If there exist Xi, . . . ,Xk G Xi and Xk+i G Xp sitc/i that Xi G Di(xo),X 2 G 
D 2 (^i), • • • G t/ien #ie event £/ + (xo) occurs. 
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Proof. • The inclusion D^y) C Cj is clear for every % G {1,...,A; + 1}. Let % G 
{2, . . . , k + 1}, y G Ci_! and z G Di(y). Then, as soon as d is large enough, 

\\z-y\\ 2 = \\z'-y'\\ 2 + \\z"-y"\\ 2 

< 2 + M //||2 + || z //||2_ 2< y " )Z "> 

d 
2 

< - + - 2d" 1 ) 2 + (d* - 2d -1 ) 2 - 2(d*_i - 3d~ 1 )(d i - 3d _1 )cos#i 

< d- + d 2 _ x - 2d i _ 1 (d i _i + din) - 2d~ l/2 d i d i ^ l + O^d" 1 ) 

< r 2 - 2d~ 1/2 d i d i ^ 1 + Oi(d" 1 ) < r 2 . 

Let now y G C and z G Di(j/). As c?i = 1 + p = r x > 2, we obtain, for d large enough : 

|| z - = || z ' - y'f + - y"f < £ + - 2d" 1 ) 2 < r 2 . 

• The second point is a simple consequence of the first point, of the fact that the sets 
DA as the sets Cj, are disjoint and of the definition of the event Q + . □ 

Note that for i G {1, . . . , k + 1}, \Di{y)\ and \D"(y")\ do not depend on the choice 
of y G Cj_i. We thus denote by and \D"\ these values. We now give asymptotic 
estimates for these values : 

Lemma 2.13. For every i G {2, . . . , k + 1}, 



1 l-D-'l 1 IDA I 

lim —In — — = lim —In — — = Infra/ 1 — a 2 ). 

d^+oo d Vd-2 d^,+oo d Vd V 

Proof. We have, by homogeneity and isotropy: 

(27) \D"\ = (fa - 2d~ 1 ) d ~ 2 - (d, - 3d~ l ) d ~ 2 ) \S\ 

where S — {x = (x\, . . . , x d ^ 2 ) £ B"(0, 1) : x\ > \\x\\ cos(#j)}. 
But S is included in the cylinder 

{(x i ) 1 < i < d . 2 eR d - 2 : Xl G [0,1], \\(x 2 ,...,x dr . 2 )\\ <sin(^)} 

and S contains the cone 

{(xi)i<i< d - 2 G M d ~ 2 : X! G [0,cos(^)], || (x 2 , . . . , x d - 2 ) \\ < x x sin(^) cos^)" 1 }. 

Therefore : 

(28) ^_ 3 cos(^)sin(ft t ) d " 3 (d-2)- 1 < \S\ < v d . 3 sin^)'*" 3 - 

From fl27|) . fl28|) . and the limits cos(#i) — > (dj_i H-G^r^d,^ 1 ^ and d, sin(0j) — >■ rj^/T 
we get 

1, 



lim -In j -^- L = Infra/ 1 - a 2 )- 

d \V d -2j v 

The lemma follows. Note that a direct calculus with spherical coordinates can also give 
the announced estimates. □ 

Everything is now in place to prove Proposition 12.91 
Proof of Proposition 12. 9L Choose p < 1 and x G W such that x" = 0. 
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• We start with a single individual, encoded by its position ( = {x} C C , and we 
build, generation by generation, its descendance : we set, for 1 < % < k, 

Ci = XiH |J Di(y)cCi, 

and for the (k + l)-th generation, we finally set 

Cfc+i = X P n [J D k+1 (y) C C fc+1 . 

■yea 

By Lemma \2.12\ if ( k +i ^ then the event Q + (x) occurs. To bound from below the 
probability that C,k+x 7^ 0, we now build a simpler process £, stochastically dominated by 
77. 

• We set a, = Ai|£>j| for i £ {1, . . . , fc} and a k +i = X p \D k+1 \ : thus, aj is the mean 
number of children of a point of the (i — l)-th generation. Let X = x be the position 
of the first individual. Let N be a Poisson random variable with parameter a\ : This 
random variable gives the number of children of Xq. 

Then we complete X in a family X = (Xq, X±, . . . , X^+i) of points in M. d , independent 
of N, in the following manner : X\ is taken uniformly in D\(Xq), then X2 is taken 
uniformly in D2(X\), and so on : The family X represents the positions of the possible 
successive descendents of x in a single branch. Let then (X^)j>i be independent copies 
of X. We will use the N first copies, one for each child of Xq. 

Let Y = (Y/)2<i<k+i,j>i be an independent family of independent and identically 
distributed random variables, such that Yf follows the Bernoulli law with parameter 
1 — exp(— aii), which is the probability that a Poisson random variable with parameter 
ai is different from 0. The random variable Y? encodes the fact that in the j-th branch, 
the member of the (i — l)-th generation has children or not in a process in which we 
forget dependencies due to geometry : Y? — 1 if and only if the member of the (i — l)-th 
generation in the j-th branch has at least one child. In other words, with the Y? , we only 
keep, for each individual from the first generation X±, at most one child, and its position 
is given by X\. We set J\ = {1, . . . , N} and, for every % 6 {2, . . . , k + 1} : 

Ji = {l<j<N:Yj = --- = Y? = 1}. 

Thus the random set Jj gives the exponents of the branches that are, among the N 
initial branches, still alive at the i-th generation in a process with no dependecies due to 
geometry. 

Until now, we did not take into account the geometrical constraints between individuals. 
For every i £ {2, . . . , k + 1} and every j > 1, we set 

Z{ = 1 if X{ £ [J D { (X{_^ and Z\ = otherwise. 

feJi-i\{j} 

When building generation i from generation i — we explore the Poisson point processes 

in the area lJ JgJ _ i By construction of the Cj, these areas are distinct for 

different generations. When Z\ = 1, we are sure that the member of the i-th generation 
in the j'-th branch does not live in an already explored zone. Thus, one can check that, 
for every % £ {2, . . . , k + 1}, the set 

£ i = {X> i :jeJ i wdZi = .-- = Zi = l} 

is stochastically dominated by Q. Thus to prove Proposition I2.9[ we now need to bound 
from below the probability that is not empty. 
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• Let T be the smallest integer j such that = • • • = Yj!+i = 1 : in other words, 
T is the smallest exponent of a branch that lives till generation k + 1. To ensure that 
^ 0, it is sufficient that T < N and that = ■ ■ ■ = Z^ +l = 1. So : 

l-P(£+(x)) < Pfe + i = 0) 



< p(#j fc+1 = o) + p({T<iV}n |J R T = o}J 

V 2<j<fc+l / 

< P(#J fe+1 = 0)+ ^ P(T < iV and = 0). 



2<i<fc+l 

For every 2 < i < k + 1, we have by construction : 

P(T < N and = 0) = P (T < N, 3j G J^a \ {T} such that Xf G A(^ti)) 

< P (T < iV and j G P_i \ {T} and G A(*ti)) 
j">i 

= ^P(Wl^_ A{T} P(Xf G A(*£-i) \Y,N)) 

= ^ p ( WWxxm) p (*< e A(xtx)) 

< P(#P_i)P (X, 1 g A(XtO) • 

Besides, as {X})" is uniformly distributed on Cf and is independent of {Xf_-^" 

(29) PKeA(4i)) = P {{Xlf e D'^XUf)) J-^ 
This leads to 

(30) 1 - P(£ + (x)) < P(# J fe+1 = 0) + ^ P(# J, 



IP" I 



i=2 



For 1 < % < k + 1 , the cardinality of Jj follows a Poisson law with parameter 

i 



i'=2 



Remember that ai = Ai|Pj| for i G {1, . . . , k} and a^+i = X p \Dk + i\. By lemma [2.13[ we 
have the following limits: 

1 . k(1 + p) 
hm — In a\ = In < 0, 

lim -lnaj = ln(/t-»/ 1 — a^) < for 2 < i < k, 

d—>-+oo a 



-1 + p, 



lim -lna fc+ i = ln(Wl-a| x — — ) < 0. 

d^+oo a V 

To see the signs of the limits, note that k < 1, that < 1 and that ff2^|) implies that 

fc+i 4 P 2 
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Consequently, we first see that 

lim -ln(?7 fe+ i) = lim -]n(an 

d— >+oo a d— »+oo a 




> with d22P; 



(31) therefore, lim P(#J fc+ i = 0) = 0. 

d— >+oo 

Similarly, for 2 < i < k + 1, we have 
lim -ln(^_i) = 

d— >+oo a 

Lemmas 12.111 and 12.131 ensure that : 




d-H-oo d VIC 

Thus, for 2 < % < k + 1, we have : 



-; In ■ „?, ■ ^ = In ^ 



Now, 



(32) for 2 < i < k, lim sup - In ( £(#</;_] 

d— i-+oo " 



(33) limsup^lnfa#J fc )|gH) < In I ^ ' " - / 




with ([21]). To end the proof, we put estimates (|5I|), (|32j) and ([33]) in (j30|). □ 

2.4.2. Several steps in the 2- dimensional oriented percolation model. We prove here Propo- 
sition [2^8] by building the supercritical 2-dimensional oriented percolation process embed- 
ded in the two-type Boolean Model. 

Proof of Proposition 12.81 We first define an oriented graph in the following manner: 
the set of sites is 

S = {(a,n)6ZxN: \a\ < n, a + n is even }; 

from any point (a,n) G S, we put an oriented edge from (a,n) to (a + l,n + 1), and 
an oriented edge from (a,n) to (a — 1,71+ 1). We denote by p c (2) G (0, 1) the critical 
parameter for Bernoulli percolation on this oriented graph - see Durrett 84 [5] . 
For any (a, n) G S, we define the following subsets of M. d 

W a<n = d" 1/2 (]fl-l,a + l[x]n-l,n[xl d - 2 ), 
W~ n = r 1/2 (]a-l,o[x]n,n + l[xR rf - 2 ), 
W+ n = r 1/2 (Ka+l[x]n,n + l[xl d - 2 ). 
Note that the (W a , n )^ n ) &s are disjoint and that W+ n U W~ +2>n C W a+hn+ i. 
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We now fix k > 1 and n G (K c p (k), 1), and for x G W a ^ we introduce the events : 

g+ ( x \ / There exist distinct Xi, . . . , Xk G Xi H W+ n and rcfc+i G Xp H W 7 ^ 1 
a ' n | such that x , Xi, . . . , x^+i is a path J ' 

g_ / s f There exist distinct ii,...,Xfc G Xi H W~ n and Xfc+i G Xp H W 7 ^ 1 

a,ri | such that x , Xi, . . . , x^+i is a path J ' 

Note that ^ ( a; ) * s exactly the event Q + {x) introduced in fl23l) . and that the other events 
are obtained from this one by symmetry and/or translation. 

Next we choose p G (p c (2), 1). With Proposition I2.9[ and by translation and symmetry 
invariance, we know that for every large enough dimension d, for every (a, n) G S, for 
every x G W a>n : 

(34) PQSinW) > P- 

We fix then a dimension d large enough to satisfy (1341) . We can now construct the random 
states, open or closed, of the edges of our oriented graph. We denote by oo a virtual site. 

Definition of the site on level 0. Almost surely, Xp H Wo,o 7^ 0- We take then some 
x(0, 0) G x P n W ,o. 

Definition of the edges between levels n and n + 1. Fix n > and assume we 
built a site x(a, n) G W 0| „ U {c>o} for every (a,n) G S*. Consider (a,n) G S 1 : 

- If x(a,n) = 00 : we decide that each of the two edges starting from (a,n) is open 
with probability p and closed with probability 1 — p, independently of everything 
else; we set z~(a,n) = z + (a,n) = 00. 

- Otherwise : 

- Edge to the left-hand side : 

- if the event Q~ n (x(a,n)) occurs : we take for z~(a,n) some point Xk+i G 
W~ n C W a -i, n +i given by the occurrence of the event, and we open the edge 
from (a, n) to (a — 1, n + 1) ; 

- otherwise : we set z~(a,n) = 00 and we close the edge from (a,n) to (a — 

1,71 + 1). 

- Edge to the right-hand side : 

- if the event Q+ n (x(a,n)) occurs : we take for z + (a,n) some point Xk+i G 
W£ n C W a+ i tn+ i given by the occurrence of the event, and we open the edge 
from (a, n) to (a + 1, n + 1) ; 

- otherwise : we set z + (a, n) = 00 and we close the edge from (a,n) to (a + 

1,71 + 1). 

For (a,n) outside S, we set z ± (a, n) = 00. 

Definition of the sites at level n+1. Fix n > and assume we determined the state 
of every edge between levels n and n + 1. Consider (a, n + 1) G S 1 : 

- If 2 + (a — 1, n) 7^ 00 : set x(a, n + 1) = z + (a — 1, n) G W 0) „ + x. 

- Otherwise : 

- if z - (a + 1, n) 7^ 00 : set x(a, n + 1) = z~(a + 1, n) G W a , n +i, 

- otherwise : set x(a,n + 1) = 00. 

Assume that there exists an open path of length n starting from the origin in this 
oriented percolation : we can check that the leftmost open path of length n starting 
from the origin gives a path in the two-type Boolean model with n alternating sequences 
of k balls with radius 1 and one ball with radius p. Thus, percolation in this oriented 

22 



percolation model implies percolation by /c-alternation in the two-type Boolean model. 
Let us check that percolation occurs indeed with positive probability. 

For every n, denote by T n the ex-field generated by the restrictions of the Poisson point 
processes xi an d Xp to the set 

d-1/2 ( R x (_oo ?r2 ) x R d ~ 2 ) . 

By definition of the events Q - remember that the (W / a ,n)(a,n.)G5 are disjoint - and by (I34I) . 
the states of the different edges between levels n and n+1 are independent conditionally to 
T n . Moreover, conditionally to J-* n , each edge between levels n and n + 1 has a probability 
at least p to be open. Therefore, the oriented percolation model we built stochastically 
dominates Bernoulli oriented percolation with parameter p. As p > p c (Z 2 ), with positive 
probability, there exists an infinite open path in the oriented percolation model we built; 
this ends the proof of Proposition 12.81 □ 



2.5. Proof of Theorem ll.4L We first prove how Propositions l2.2l and 12.81 give ffTTj) when 
a = 1, b > 1 and a — (3 — 1, and then we see how we can deduce the general case by 
scaling and coupling. 

When a = 1, b > 1 and a = (3 = 1. Set p = b. In this case, v = 5i +S P , so u d = 5i~\ — -:S P . 



Note then that the two-type Boolean model E introduced in Subsection 12.11 and whose 
intensities depend on k e (0, 1) coincides with the Boolean model directed by the measure 



K d 



v d 2 d 



I'd 



as defined in the introduction. 

If k < n c p then, by Proposition [221 there is no percolation for d large enough. Therefore, 
for any such k and for any large enough d we have: 



X c d (y d ) > — A and then \ c d {u d ) = X d (u d )v d 2 d J r d v d {dr) > 1k a . 
Letting d goes to +oo and then k goes to n c p , we then obtain 

(35) liminfiln(A>,))>ln(<). 

As k c p < 1 by Lemma |2.1[ choose now k such that k c < k < 1. Then, there exists 
k > 1 such that K c p (k) < k. Therefore, by Proposition I2.8[ there is percolation for d large 
enough in S; by coupling, this remains true for any k > K c p . Therefore, for any k > K c 
and for any large enough d we have, as before: 

A^)<^4 and then \ d (is d )<2n d . 
v d 2 a 

Letting d goes to +oo and then k goes to k c , we then obtain 

(36) limsup^-ln(A^)) < In (^) . 

d— ¥+OD U> 



Bringing f )35|) and fl36|) together, we get (JTTJ when a — 1, b — p > 1 and a — (3 — 1. 
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When b > a > and a = j3 — 1. Set p = b/a. Here, z/ = 5 a + set /x = 5% + 5 p . With 
the notation of the introduction, 

By the scaling relations (J2J) and ([3]), we obtain 

The result when 6 > a > and a = (3 = 1 follows then from the previous case. 

When b > a > and a, j3 > 0. Here z/ = a<5 a + Set /i = 5 a + 5&, m = min(a, (3) and 
M = max(a, /?). Then mp d < < M/i^ and so 



y r d dp d (r) < J r d dv d {r) < M J r d dp d (r) 



jf X d(Vd) = X&Mm) < \ c d {v d ) < X d (mfi d ) = ^X C M. 

The two previous inequalities give: 

m ~ ~ M ~ 

M m 
and the theorem follows from the previous case. □ 

3. Proof of Theorem 11.21 

Theorem 11.21 follows from Theorem 11.41 by coupling and scaling. By assumption, /i is a 
measure on (0, +oo) whose support is not a singleton. We can therefore choose b' > a' > 
in the support, set p = b'/a' and then take a small enough e > such that 

a'(l + e) < b'(l - e), p([a'(l - e), a'(l + e)] > 0, 
[i([b'(l - e), b'(l + e)] > 0, (1 + e)(l - e)^^ < 1. 

Set a = a'(l-e),b = b'(l-e) and r = (l + e)(l- £ )- 1 > 1. We have 

ar < b, p([a,ar]) > 0, p([b,br]) > and rn c p < 1. 

Set v = /i([a, ar})5 a + p([b, br])Sb and S = [a, ar] U [b, br]. For all d > 1 we have 

T~ d u d ({a}) = p([a,aT})(aT)~ d < / r~ d p(dr) = fi d ([a, ar}) 

J [a,ar\ 

and, similarly, T~ d u d ({b}) < p d ([b, br]). By coupling, this implies that X d (lsPd) < 
X d {r~ d u d ) 1 and then that 

X c d (p d ) < X d (l s p d ) < X d (T~ d u d ) = r d X d (u d ). 

But Xftfid) = Ad(yUd)2%/i((0, +oo)) and, similarly, X c d {v d ) = X c d (v)2 d v d is((0, +oo)), which 
leads to 

W„ n <r ^ M(0»+oo)) Tc^, ^ 
ArfM " T ,((0,+oo)) Arf(l/d) - 

But by Theorem 11.41 we have 

lim -In (x c d (u d )] = ln(/c„), and then limsup^ln (ASf/^)) < ln(r«„) < 0, 

d^+oo a \ J p d^+oo d \ J p 

which ends ce proof. □ 
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Note that as a byproduct of the proof, we obtain the following upper bound : 
limsup^ln CxiQid)) < inf ln(^ /a ). 

>+oo U ^ / 0<a<b<+oo, a,oebupp(/i) 

An interesting but certainly tough result would be to obtain a real limit, and there is no 
reason why this limit should coincide with this upper bound. 

References 

[1] Ajit Balram and Deepak Dhar. Scaling relation for determining the critical threshold for continuum 
percolation of overlapping discs of two sizes. Pramana, 74:109-114,2010. 10.1007/sl2043-010-0012-0. 

[2] R. Consiglio, D. R. Baker, G. Paul, and H. E. Stanley. Continuum percolation thresholds for mixtures 
of spheres of different sizes. Physica A: Statistical Mechanics and its Applications, 319:49 - 55, 2003. 

[3] Deepak Dhar. On the critical density for continuum percolation of spheres of variable radii. Physica 
A: Statistical and Theoretical Physics, 242(3-4):341 - 346, 1997. 

[4] Deepak Dhar and Mohan K. Phani. Continuum percolation with discs having a distribution of radii. 
J. Phys. A, 17:L645-L649, 1984. 

[5] Richard Durrett. Oriented percolation in two dimensions. Ann. Probab., 12(4):999-1040, 1984. 

[6] Jean-Baptiste Gouere. Subcritical regimes in the Poisson Boolean model of continuum percolation. 
Ann. Probab., 36(4):1209-1220, 2008. 

[7] Jean-Baptiste Gouere. Subcritical regimes in some models of continuum percolation. Ann. Appl. 
Probab., 19(4):1292-1318, 2009. 

[8] Jean-Baptiste Gouere. Percolation in a multiscale boolean model. Preprint, 2010. 

[9] Janos Kertesz and Tamas Vicsek. Monte carlo renormalization group study of the percolation prob- 
lem of discs with a distribution of radii. Zeitschrift fiir Physik B Condensed Matter, 45:345-350, 
1982. 10.1007/BF01321871. 

[10] Harry Kesten. Asymptotics in high dimensions for percolation. In Disorder in physical systems, 
Oxford Sci. Publ., pages 219-240. Oxford Univ. Press, New York, 1990. 

[11] Ronald Meester and Rahul Roy. Continuum percolation, volume 119 of Cambridge Tracts in Math- 
ematics. Cambridge University Press, Cambridge, 1996. 

[12] Ronald Mecstcr, Rahul Roy, and Anish Sarkar. Nonunivcrsality and continuity of the critical covered 
volume fraction in continuum percolation. J. Statist. Phys., 75(1-2):123-134, 1994. 

[13] M. V. Menshikov, S. Yu. Popov, and M. Vachkovskaia. On the connectivity properties of the com- 
plementary set in fractal percolation models. Probab. Theory Related Fields, 119(2):176-186, 2001. 

[14] Mathew D. Penrose. Continuum percolation and Euclidean minimal spanning trees in high dimen- 
sions. Ann. Appl. Probab., 6(2):528-544, 1996. 

[15] John A. Quintanilla and Robert M. Ziff. Asymmetry in the percolation thresholds of fully penetrable 
disks with two different radii. Phys. Rev. E, 76(5):051115, Nov 2007. 

Laboratoire de Mathematiques, Applications et Physique Mathematique d'Orleans 
UMR 6628, Universite d'Orleans, B.P. 6759, 45067 Orleans Cedex 2 France 
E-mail address: Jean-Baptiste.Gouere@univ-orleans.fr 

Institut Elie Cartan Nancy (mathematiques), Universite Henri Poincare Nancy 1, 
Campus Scientifique, BP 239, 54506 Vandoeuvre-les-Nancy Cedex France 
E-mail address: Regine.Marcliand@iecn.u-nancy.fr 



25 



